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ABSTRACT 

Several numerical simulations of the transient flow of helium in an expansion tube are 
presented in an effort to identify some of the basic mechanisms which cause the noisy test 
flows seen in experiments. The calculations were performed with an axisymmetric Navier- 
Stokes code based on a finite-volume formulation and upwinding techniques. Although 
laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed 
some of the important features seen in the experiments. In particular, the discontinuity 
in tube diameter at the primary diaphragm station introduced a transverse perturbation 
to the expanding driver gas and this perturbation was seen to propagate into the test gas 
under some flow conditions. The disturbances seen in the test flow can be characterized as 
either small-amplitude, low-frequency noise possibly introduced during shock compression 
or large-amplitude, high-frequency noise associated with the passage of the reflected head of 
the unsteady expansion. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science and 
Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. The author’s present address: 
Department of Mechanical Engineering, The University of Queensland, St Lucia, Qld 4072, AUSTRALIA. 
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Nomenclature, Units 


A 

a 

c P 

c v 

E 

e 

e 

F 

G 

h 

k 

h 

P 

Pr 

Q 

9 

R 

Re 

r 

T 

t 

U_ 

u 

u 

v 

X 

y 


(x,y)-plane cell area 

local speed of sound, m/s 

specific heat at constant pressure, J /kg 

specific heat at constant volume, J/kg 

total energy (internal + kinetic), J/kg 

specific internal energy, J/kg 

unit vector 

algebraic vector of x-component fluxes 
algebraic vector of y-component fluxes 
specific enthalpy, J/kg 
coefficient of thermal conductivity 
unit normal vector 
pressure, Pa 

Prandtl number, (C p g/k) 
algebraic vector of source terms 
heat flux 

gas constant, J/kg/K 
Reynolds number 
radial coordinate, m 
temperature, K 
time, s 

algebraic vector of conserved quantities 
Riemann invariant, m/s 
x-component of velocity, m/s 
y-component of velocity, m/s 
x-coordinate, m 
y-coordinate (radial), m 


P 

P 

7 

A 

r 

P 

v 

ft 

ft' 


: parameter in Mirels’ boundary layer theory 
: density, kg/m 3 
: ratio of specific heats 
: second coefficient of viscosity 
: shear stress, Pa 
: coefficient of viscosity, Pa.s 
: frequency, s -1 
: cell volume, m 3 

: volume per radian for the axisymmetric cell 


< • > : cell-averaged value 
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Subscripts 

0 : conditions immediately behind the shock 

1 ... 20 : gas state as shown on Fig. 1 

e : conditions just outside the boundary layer 

i : interaction of expansion with driver-gas/test-gas interface 

v : viscous contribution 

w, wall : wall value 

x, y : cartesian components 
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1 Introduction 


The expansion tube is a member of the family of pulse-type aerodynamic facilities designed 
to provide high enthalpy test gas for short durations. It has a fundamental advantage 
over the reflected-shock tunnel because its operating cycle does not involve the stagnation 
of the heated test gas. Hence, it can provide a final test flow with low ionization and 
low dissociation. Initial theoretical studies of the expansion tube [1] indicated that such 
a machine could provide a wide range of test flows by simply altering the initial filling 
pressures. However, operational experience (see e.g. [2, 3]) indicated that the test flow was 
often contaminated by large-amplitude disturbances. Here, some numerical simulations of 
an entire (but idealized) expansion tube are presented in order to provide some insight into 
the mechanisms causing the disturbances in the test flow. The specific facility configuration 
and the flow conditions considered here are based on the experimental perfect-gas study 
undertaken by Shinn and Miller [3], 

1.1 Facility Operation 

An expansion tube is essentially a single tube divided into three sections by “primary” and 
“secondary” diaphragms as shown in lower part of Fig. 1 (which is reproduced from [3]). The 
strong primary diaphragm separates the driver tube and the intermediate (or shock) tube 
while the light secondary diaphragm separates the intermediate tube from the acceleration 
tube. The intermediate and acceleration tubes have the same diameter. Initially, the driver 
tube contains a “driver gas” at high pressure (state 4) while the intermediate tube contains 
the test gas (state 1) and the acceleration tube contains a low pressure “acceleration gas” 
(state 10). The aerodynamic model to be tested is located near the downstream end of the 
acceleration tube. 

The operation of the facility, starting with the rupture of the primary diaphragm at t = 0, 
is shown schematically in the top part of Fig. 1. The numbering of the flow states is the same 
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as that introduced by Trimpi [1]. High pressure driver gas expands into the intermediate 
tube (state 3) and shock-compresses the test gas (state 2). Simultaneously an expansion fan 
travels upstream into the driver tube. On reaching the end of the intermediate tube, the 
shock ruptures the secondary diaphragm (ideally without producing disturbances) and the 
leading portion of the shock-compressed test gas is allowed to expand into the acceleration 
tube to reach the test flow conditions (state 20). The expanded test gas will have a high 
velocity (typically 6 — 8 kmjs ) and a relatively low temperature (several hundred degrees 
Kelvin). 

As shown in Fig. 2, test time commences after the passage of the test-gas/acceleration- 
gas interface and ideally finishes with the arrival of the downstream end of the unsteady 
expansion (and a corresponding increase in Pitot pressure). Not all of the shock-compressed 
test gas is expanded to the test flow state (20) and so test times are short, typically measured 
in tens to hundreds of microseconds. Unfortunately, ideal operation of the facility is seldom 
realized as test flows are typically very noisy and may have peak-to-peak variations of 50% 
in Pitot pressure. Note that, depending on the length of the tube sections and the initial gas 
states, the arrival times for the various waves at the test-section may change. Of particular 
importance is the reflected-head of the unsteady expansion (of the test gas). This wave is 
generated when the upstream-head of the unsteady expansion reaches the driver-gas/test- 
gas interface. A partial reflection will occur and, for some test conditions, this wave (which 
travels along a. u + a characteristic) will arrive at the end of the acceleration tube before the 
unsteady expansion itself. 

1.2 Previous Studies 

The cause of the test flow disturbances has been the focus of several studies. Shinn and Miller 
[3] made experimental measurements in the NASA Langley facility using helium as the test 
gas, driver gas and acceleration gas in order to eliminate chemical effects. They observed that 
most of the operating conditions resulted in test flows which were very noisy but identified 
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an optimum condition (Pi = 3.45 kPa , Pi 0 = 16 Pa) with approximately 300 fis test time 
(with only moderately small disturbances) before the onset of large-amplitude fluctuations 
in the Pitot pressure. A large dip in the Pitot pressures was observed for values of P 10 
significantly larger than 20 Pa and this was attributed to boundary layer transition. They 
also identified the importance of the diaphragm dynamics on the quality of the test flow, 
especially at low initial pressures. However, they did not identify the mechanism causing 
the early disturbances in the test gas. 

More recently, equivalent experiments [4] using both helium and argon as the working gas 
were performed at the University of Queensland to check the Langley results [3] and see if the 
useful range of test conditions could be expanded. In conjunction with these experiments, a 
numerical study (based on quasi-one-dimensional modelling) was undertaken [5] to see if the 
test gas disturbances could be caused by “blobs” of driver gas being accelerated through the 
expanding test gas. It was concluded that, even if the blobs existed, they would not arrive 
early enough to cause the test-flow disturbances. However, the reflection of the unsteady 
expansion off the driver-gas/test-gas interface could (under some conditions) arrive early 
enough to introduce the noise. Test time estimates based on this mechanism are given in [6] 
and are consistent with the present numerical simulations. 

Pauli [4, 7, 8] has pursued the possibility of acoustic disturbances being generated near 
the primary diaphragm station and then being propagated through the unsteady expansion 
of the driver gas. A family of transverse waves can exist and, when propagated through the 
expansion, will undergo a frequency shift such that a single dominant frequency will emerge. 
Transmission of these waves across the interface and into the test gas is assumed to be 
determined by the acoustic impedance mismatch at the interface. For situations where the 
sound speed in the test gas is much larger than that in the expanded driver gas, attenuation 
of the waves in the test gas should occur. Otherwise, a significant fraction of the noise 
will be transmitted. This acoustic model is seen to be a reasonable approximation to the 
simulations presented here. 
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2 The Numerical Method 


The computations reported here were performed with a finite-volume upwind code based 
on the full Navier-Stokes equations. The code is described briefly here but further detail is 
available in [9]. 

2.1 Governing Equations 

For an axisymmetric flow (with y as the radial coordinate), the finite-volume formulation of 
the Navier-Stokes equations may be expressed as 

+ ^ !s (yF - yK) dy ~b L {yG ~ yG " ] dx = Q ' (1) 


where < • > indicates a volume average, f Y is the volume per radian and S is a contour in 
the (x,y)-plane around the volume. The [/, F and G vectors are 
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Treating the viscous contributions in the form yr avoids any difficulties with the geometry 
singularity at y = 0. The effective source term is 

, „ v . / du dv v\ /cA 

, where Tgg = 2 p — \- A [ — — f — I — • (8) 

V \dx dy yj 

To augment these equations, the equation of state for a calorically perfect gas is used 

P = p( 7-l)e , ( 9 ) 

with 7 = 1.667 and R = 2077 J/kg/K for helium. The corresponding specific heats are 
C v = 3114 J/kg/K and C p = 5191 J/kg/K. The first coefficient of viscosity is evaluated 
as 

p = 5.023 x 10 -7 T 0647 Pa.s , (10) 

and the second coefficient is obtained from Stokes’ hypothesis A = — \p. Also, a constant 
Prandtl number of 0.67 is assumed. 

2.2 Numerical Implementation 

The flow domain in the (x, y)-plane is discretized as a structured mesh of quadrilateral cells 
with flow properties stored at the cell centres. At each time step, the inviscid-flux vectors 
(2) are evaluated by first applying a generalized MUSCL interpolation scheme [10] to obtain 
“left” and “right” states at the midpoints of the cell interfaces. A locally one-dimensional 
(approximate) Riemann solver is then applied to obtain the interface flow properties during 
the time step. The spatial derivatives used in the viscous flux vectors (3)-(7) are obtained 
at the cell vertices by applying the divergence theorem. The source term (8) is evaluated at 
the cell centres. The line integrals in (1) are then evaluated using the midpoint rule and the 
solution advanced in time. 

Although a predictor-corrector time-stepping scheme was available in the code, the solu- 
tions shown in the following sections have used Euler time stepping in order to reduce the 


Q' = 


0 

0 

{ P - Tgg)A/Q' 
0 
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required cpu time for each simulation. Several simulations of the initial shock compression of 
the test gas were run to investigate the effect of grid resolution and choice of time-stepping 
scheme [11]. Except for the fine-scale features, there were few differences in the solutions for 
both Euler time-stepping and predictor-corrector time-stepping on a (coarse) 1604 x 20 grid. 
The differences between these coarse-grid solutions and the solution on the (fine) 2406 x 30 
grid were also small. The limit on the number of cells was set by the memory available on 
the local workstations (64 Mbytes) and the desire to fit the jobs within the 8 Mword queues 
on a Cray-YMP located at the NASA Langley Research Centre, USA. 
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3 Simulation Results 


Figure 3 shows the flow domain considered in the simulations. The flow is considered to 
be axially symmetric with the flow domain extending radially from the axis to the tube 
wall and axially from the closed end of the driver tube to the end of the acceleration tube 
where it enters the test section of the facility. Dimensions are shown in the figure. Using 
the program in [11], the grid was generated in three pieces, one for each section of the 
facility. The number of cells and the clustering toward the tube wall remained constant 
along the length of the domain and the cell size remained constant in the axial direction. 
The radial distance from the outermost cell-centres to the intermediate- and acceleration- 
tube walls is 3.87 X 10 -4 m. Although the use of an axisymmetric geometry minimizes 
the computational resources required for the simulation, it also precludes genuinely three- 
dimensional disturbances which may constitute a significant portion of the experimentally 
observed noise. 

While the change in tube radius at the primary- diaphragm station is included, neither 
the diaphragm nor the square transition section downstream of this point is modelled. In- 
stead, the simulations start with high-pressure reservoir conditions upstream of the primary 
diaphragm station and intermediate-tube initial conditions downstream of the same station 
(and into the acceleration tube). The simulations are allowed to proceed until the primary 
shock reaches the secondary-diaphragm station and, at that point in time, the conditions 
downstream of the secondary diaphragm are reset to the acceleration-tube initial conditions. 
Thus, the bursting of the diaphragms is considered to be ideal in the following calculations. 
This is in contrast to the experimental observation that the flow quality is very sensitive to 
rupture behaviour of the diaphragms [2, 3] but, as will be seen in the results, the key features 
governing the quality of the test flow are captured. 
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3.1 Initial State 


All of the simulations used the nominal driver tube conditions given in [3]. These are 

P 4 = 33 MPa, u 4 = 0, T 4 = 330 K, p 4 = 48.15 kg/m 3 , e 4 = 1.028 MJ/kg. 

As discussed in [3] the nominal driver-tube temperature may be lower than the actual gas 
temperature at primary-diaphragm rupture where a maximum value of T 4 = 390 K was 
estimated. The test gas and acceleration initial conditions are 

T i0 = T\ = 300 K, eio = e 4 = 0.934 MJ/kg, uj 0 = itj = 0, 

with pressures and densities as specified in Table 1. 

For cases la and lb, the conditions were chosen to approximate the optimum conditions 
identified in the experiment [3] and two different grids were used to check for sensitivity to 
grid resolution. Cases either side of this optimum condition were then considered. Case 2 has 
a higher intermediate-tube pressure where large-amplitude noise was observed at an earlier 
time in the experiment. Case 3, with a low intermediate-tube pressure, also was observed 
to have a noisy test flow early in the test time. Finally, case 4 considers the condition with 
high acceleration-tube pressure where (for the Langley expansion tube) a large dip in Pitot 
pressure was observed shortly after the arrival of test flow [3]. 

3.2 Case 1: P\ = 3.45 kPa , P 10 = 16 Pa 

This condition was chosen to be the starting point for the study because it was identified 
as the best operating condition in the experiment [3]. The observed test time was approx- 
imately 300 ps before large-amplitude noise contaminated the test flow. During the test 
time, moderately small disturbances in Pitot pressure were observed. Two calculations were 
performed and details of the computational resources are shown in Table 2. The global time 
step was limited so that the maximum Courant-Friedrichs-Lewy (CFL) number in any cell 


10 



was approximately 0.5. This maximum value occurred in the flow region just behind the 
shock. 

Case la used a “coarse” grid of 1604 cells axially by 20 cells radially while case lb used 
a “fine” grid of 2406 x 30 cells. (Refer to [11] for the coarse grid results.) The computed 
results, including the timing of events such as the passage of the shock and contact surfaces, 
were largely the same. There were, however, some differences in the details of the contact 
surfaces and the noise introduced to the test flow. In particular, the coarse grid attenuated 
some of the noise and limited the highest frequency seen in the simulation. 

Figure 4 summarizes the calculation for case lb as a wave diagram. The data points 
plotted were obtained from the individual plots of the density field shown later. Note the 
small size of the slug of compressed acceleration gas (state 20). The position of the long 
dashed line, denoting the downstream-tail of the unsteady expansion of the test gas, was 
estimated by its arrival at x = 24 m while the short-dashed line represents the propagation 
of the leading edge of the reflected expansion and the asterisk denotes the arrival of large- 
amplitude noise at the end of the acceleration tube (a; ~ 24 m). This reflected expansion is 
the result of the interaction between the upstream-head of the unsteady expansion and the 
test-gas/ driver-gas interface. The test time (state 5) for this case starts after the passage of 
the test-gas/acceleration-gas interface and stops with the arrival of the downstream tail of 
unsteady expansion. 

3.2.1 Shock Compression of the Test Gas 

Density contours for the shock compression of the test gas are shown in Fig. 5. For each 
frame, 81 contours have been plotted over the range —6.0 < log 10 /> < 2.0. Although the 
shock has remained essentially planar, the driver-gas/test-gas interface has become quite 
distorted by t = 2.0 ms. Boundary layer effects are present (as can be seen by the contact 
surface being left behind along the tube wall) but they are not well resolved in this part of 
the calculation. 


11 



Of special interest is the train of transverse waves following the interface. These transverse 
waves show up clearly in the contours of the radial velocity as shown in Fig. 6 where it 
appears that the disturbances are introduced at the primary diaphragm station. Note that 
the axial scale is highly compressed and the actual Mach angles are much shallower than 
those that appear in the figure. Significant noise is present in the shocked test gas. 

Figure 7 shows the axial variation of flow properties at t — 1.0 ttis together with the 
estimated properties for a purely one-dimensional situation. Except for some noise, the 
levels are in general agreement. Although the axisymmetric simulation exhibits a higher 
shock speed than the one-dimensional case, tables 3 and 4 show that the computed shock 
speed is still lower than that measured in the experiment. The large experimental value was 
discussed in [3] and may be partly due to a low estimate for the initial driver-gas temperature. 
Also, note the variation in shock speed over this phase of the facility operation (see table 4). 

3.2.2 Unsteady Expansion of the Test Gas 

The unsteady expansion of the test gas is shown in Fig. 8. Viscous effects are more important 
in the acceleration tube and a very thick boundary layer is established behind the shock and 
test-gas/acceleration-gas interface. The separation distance of the shock and contact remains 
essentially constant from t = 2.5 ms to t = 4.0 ms because a balance has been reached 
between the mass flow of acceleration gas entering through the shock and the mass flow 
bleeding into the boundary layer [12]. The hot acceleration gas passing into the boundary 
layer tends to accumulate just upstream of the test-gas/acceleration-gas interface and form 
a large bulge in the boundary layer (at i ~ 23 m, t = 4.5 ms). It is not until the last frame 
(t = 5.0 ms) that any significant level of noise is evident in the expanding test gas. 

Figure 9 shows the axial variation of flow properties at t = 4.4 ms. There are very obvious 
disturbances in the region of the expanded driver gas (state 3) and, on close inspection of the 
pressure distribution, some small-amplitude disturbances can be seen near the downstream 
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end of the unsteady expansion of the test gas (at x ~ 22 m). Pauli’s acoustic wave model [7] 
suggests that the broadband disturbances, introduced upstream of an unsteady expansion, 
will tend to focus to a single frequency near the downstream tail of the expansion. For 
the expanding driver gas, Pauli estimates a “focus” frequency vp — 11.9 kHz while the 
perturbations seen in the present simulation have an approximate frequency v = 10.2 kHz. 
These figures are in reasonable agreement. 

The very high temperature of the shocked acceleration gas enhances the diffusion near 
the test-gas/acceleration-gas interface. Figure 10 shows details of the shock and test- 
gas/acceleration-gas interface region. The bulge seen very clearly at x ~ 23 m, t = 4.5 ms 
in Fig. 8 is no longer as obvious in this “true-shape” plot. Plotting the velocity field in a 
frame of reference which is stationary with respect to the shock and contact surface clearly 
shows the bleeding of shocked acceleration gas out of region 20 and into the boundary layer. 
Estimates of the separation distance between the shock and contact surface (i.e. the length 
of region 20) can be made using the theory presented in [12]. Using the expression 

Lm = 7p{ 

together with the parameter estimates 

/? = 2.12, u w = U a , io = 5900 m/s, u ei0 = 1200 m/s, 

Pe,o = 1.0 x 10 -4 kg/m 3 , p w p = 2.484 x 10~ 4 kg/m 3 , p Wt o = 2.012 x 10 -5 Pa.s , 

a separation length of 0.20 m is obtained. This is in reasonable agreement with the 0.22 m 
length observed in Fig. 10 but it should be noted that the value of L 20 obtained from eqn. 
(11) is very sensitive to the estimated parameters. 

3.2.3 Histories of the Test Flow Properties 

Histories of the centreline test-flow properties near the exit of the acceleration tube are 
shown in Fig. 11. Features that can be identified include the passage of the shock and 


Pe y 0 \ ^e,G ^e t 0 Pw } 0 

Pw , 0 / Uw ^e,0 Pw , 0 
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the test-gas/acceleration-gas interface at times t = 4.45 ms and t = 4.48 ms respectively. 
These events are followed by a gradual variation in properties possibly due to a constriction 
of the effective tube diameter associated with the passage of the bulge (identified in Fig. 
8), The downstream tail of the unsteady expansion arrives at t = 4.63 ms when average 
pressure and density begin to rise while Mach number and velocity decline. It is in the 
period 4.6 < t < 5.0 that a “small-amplitude” fluctuation is evident. The period of this 
perturbation is approximately 40 fis, which corresponds closely to the period observed in the 
Shinn and Miller [3] experiment and the more recent calibration measurements [13]. 

The fluctuations increase in both the amplitude and the frequency with the arrival of 
the reflected head of the unsteady expansion at t — 5.04 ms. This wave originates with the 
rupture of the secondary diaphragm and travels along a u 2 — «2 characteristic until it strikes 
the driver-gas/test-gas interface. There is both transmission into the driver gas and partial 
reflection downstream along a u + a characteristic thus admitting more disturbances into the 
expanding test gas. The arrival time for this wave (at the exit plane of the acceleration tube) 
may be estimated as described in [6]. The process starts at x' = 0, t' = 0 with the rupture 
of the secondary diaphragm. The upstream head of the unsteady expansion intersects the 
driver-gas/test-gas interface at t' = A t{, x' = ( u 2 — a 2 ) Afj. Refer to tables 5 and 6 for 
particular values. The propagation time of the wave through the unsteady expansion may 
be obtained by integration along a u + a characteristic for which the Riemann invariant is 
[14] 

2—2 

u 5 4 -a 5 = U = u + -a . (12) 

7-1 7-1 

Within the centred expansion we have 


V 


— u — a , 


( 13 ) 


and this may be substituted into equation (12) to obtain 


u 


(14) 


7-1 
7 + 1 



2 x'\ 


The integration may be terminated at the end of the expansion (a = a 8 ) or at the end of the 
acceleration tube, whichever is sooner. If the reflected wave reaches the downstream end of 
the expansion before it reaches the end of the tube, it continues to propagate into the test 
flow region (20) with speed 115 + 05. The time of arrival as computed with this approach 
closely matches the onset of large-amplitude noise in the exit flow histories and is indicated 
by an asterisk in the wave diagram (Fig. 4). 


A comparison of computed test flow properties with the experimental measurements of 
Shinn and Miller [3] is shown in Fig. 12. The values are taken from the Pitot histories 150 fis 
after the passage of the shock. The overall level of the computed wall pressure is a little low, 
possibly because of the low estimate for the initial temperature in the driver tube. Also, 
the computed wall pressure is smoother than the experimental wall pressure. This may be 
partly due to numerical dissipation and partly due to noisy measurement techniques in the 
experiment. The arrival of the downstream tail of the unsteady expansion can be seen by 
the rise in pressure 0.2 ms after the passage of the shock in both the simulation and the 
experiment. 


The other frames in Fig. 12 show Pitot pressures. The computed Pitot pressure at the 
centreline is again a little lower than the experimental value and is also less noisy for the 
first 0.4 ms after shock arrival, possibly due to excessive diffusion in the numerical scheme. 
However, there is good agreement in that two types of disturbances are evident: a small- 
amplitude noise before the arrival of the reflected head of the unsteady expansion and a 
large-amplitude noise after the passage of the reflected expansion. This is one of the main 
results of this study. 


Other features evident in the Pitot histories include the large dip at t = 5.0 ms, r = 
52.5 mm and the consistently low values at r = 71 mm. The large dip is associated with 
the passage of the bulge in the boundary layer while the r = 71 mm position remains within 


15 


the boundary layer after the passage of the test-gas/acceleration-gas interface. 

3.3 Case 2: Pi = 20.62 kPa , Pio = 16 Pa 

Having simulated the arrival of noise in the test flow at F x = 3.45 kPa, it was decided to 
try a higher pressure of Pi = 20.62 kPa where the experimental measurements indicated 
that large-amplitude disturbances appear very early in the test flow. Fig. 13 shows the wave 
diagram for this case. Note that the downstream tail of the unsteady expansion arrives at 
the tube exit shortly after the test flow begins and that the reflected head of the unsteady 
expansion is expected to arrive much later. Because of the higher initial intermediate-tube 
pressure, the shock has a lower speed than in case 1. 

The first frame of Fig. 14 shows the state of the density field just before the rupture of 
the secondary diaphragm. Transverse waves are still present in the expanded driver gas but 
they seem to be weaker. Also, there is some noise evident in the compressed test gas. 

The unsteady expansion of the test gas is shown in sequential frames of Fig. 14. The 
process appears similar to that in case 1 except that there is very little noise evident in the 
expanded test gas (even at t — 5.0 ms). This lack of noise is confirmed in Fig. 15. Although 
perturbations have been introduced at the primary-diaphragm station, the weaker expansion 
of the driver gas has resulted in a weaker focusing of the noise at the downstream end of the 
expansion (when compared to that in case 1). 

Fig. 16 compares the computed pressures at the exit of the acceleration tube with 
experimental measurements. The arrival of the downstream-tail of the unsteady expansion 
can be clearly seen at t ~ 4.8 ms while the reflected head of the unsteady expansion is not 
expected to arrive at the tube exit until t ~ 5.84 ms (see table 6). The lack of noise in 
the laminar solution suggests that the disturbances in the experimental test flow are caused 
by boundary layer transition or turbulence (which is not modelled in these simulations) 
rather than propagation of noise through the driver-gas/test-gas interface. The omission of 


16 


turbulence is a major limitation of this study as the experimental flow is turbulent over most 
of the high-pressure operating range of the facility. 

The average value of the computed Pitot profile (taken 150 fis after passage of the shock) 
is again slightly lower than the experimental value, possibly because of the low initial value 
for the driver gas temperature. The very low Pitot pressure values computed for r/r w ~ 0.7 
are the result of taking the value of Pitot pressure in the hot acceleration gas following 
behind the interface. 

3.4 Case 3: Pi = 1.74 kPa , Pio = 16 Pa 

Given the failure of the P\ = 20.62 kPa case to exhibit large disturbances in the simulated 
test flow, a case with lower intermediate-tube pressure than case 1 was attempted. Figures 
17 - 19 show the results of a simulation with Pi — 1.74 kPa, Pio = 16 Pa. 

Table 4 shows that the shock speed during the initial compression of the test gas is 
slightly higher than for case 1 but that the final shock speed during the unsteady expansion 
process is slightly lower. The shock compression of the test gas appears to be noisier than for 
case 1 and the centreline histories indicate that small-amplitude disturbances appear very 
early in the test flow. This seems to be consistent with the experimental Pitot pressure 
results shown in Fig. 19. The arrival of the reflected head of the unsteady expansion 
(f ~ 5.1 ms) is marked with a double asterisk in Fig. 17 and can be seen as a decrease in 
the growth rate of Pitot pressure in Fig. 19. In the simulation, the arrival of this wave is 
accompanied by larger-amplitude and higher-frequency disturbances which are not obvious 
in the experimental measurements. Also, the average value of the computed Pitot pressure 
across the exit plane is significantly lower than the experimental value. 

3.5 Case 4: Pi = 3.45 kPa , P 10 = 52.6 Pa 

For large values of acceleration-tube filling pressures, Shinn and Miller [3] also observed a 
large dip in the Pitot pressure shortly after the start of the test flow. The phenomenon 
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was attributed to boundary layer transition. This particular case simulates such an op- 
erating condition in order to see if the bulge in the boundary layer (following the test- 
gas/acceleration-gas contact surface) contributed to the large dip. 

The results for the laminar-flow simulation with P\ = 3.45 kPa , P I0 = 52.6 Pa are shown 
in figures 20 - 22. The shock compression phase is the same as for case la but, with the larger 
mass of acceleration gas, the separation distance between the shock and the contact surface 
is increased and equilibrium is attained by t = 2.8 ms. The histories of the wall pressure 
shown in Fig. 22 indicate that small-amplitude noise is present in both the acceleration gas 
and test gas up until the arrival of the reflected head of the unsteady expansion at t ~ 4.9 ms. 
This event is denoted on the wave diagram (Fig. 20) by a double asterisk The arrival time, 
estimated with equations (12) - (14) is denoted by a single asterisk. 

Although the large-amplitude noise after the passage of the reflected expansion seems to 
be consistent with the experimental measurements, there is no evidence of the large Pitot- 
pressure dip in the laminar simulation. This is consistent with the hypothesis [3] that the 
dip is caused by boundary layer transition. Note that this case was computed on a coarse 
mesh which may result in the attenuation of the high-frequency components of the noise. 
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4 Concluding Remarks 

From the computational point of view, it is possible to perform reliable simulations for 
laminar axisymmetric flow of a perfect gas in an expansion tube over a limited range of 
operating conditions. Although viscous effects are adequately resolved in the acceleration 
tube, they were not well resolved in either the intermediate tube or the driver tube. This 
is not considered to be a major problem as inviscid multidimensional flow features seem to 
dominate the performance of the driver and intermediate sections of the facility. 

The numerical technique is capable of capturing strong shocks and contact surfaces 
with large temperature jumps, however the contact surfaces tend to be more (numerically) 
smeared than the shocks. The result is that the estimates of wave propagation times and 
test flow properties appear to be in reasonable agreement with experimental measurements. 
However, the issue of reliably capturing the high-frequency noise is not fully settled and fur- 
ther computational experiments using alternative higher-order schemes (such as ENO [15]) 
are required. 

A laminar perfect-gas flow was deliberately selected to keep this study manageable, how- 
ever, there is an immediate need to perform simulations with nonideal thermal and chemical 
effects with turbulence in three-dimensions. Diaphragm dynamics is another important item 
which needs to be considered for future study. In the near future, it is unlikely that all of 
these features could be included but the modelling of turbulence and the diaphragm dy- 
namics are key items required to simulate the expansion tube over a larger range of useful 
operating conditions. 

With respect to facility design and operation, the simulations show some important 
features. Firstly, the discontinuity at the primary diaphragm station is sufficient to produce 
disturbances which may then propagate into the test flow. Secondly, the simulations clearly 
identify two categories of noise in the test flow. One is a small-amplitude, low-frequency 
disturbance which may be seen before the arrival of the reflected head of the unsteady 
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expansion while the other is a large-amplitude, high-frequency noise associated with the 
passage of the reflected expansion. 

This distinction explains the experimental observation [4, 5] that increasing the acceleration- 
tube pressure _P 10 while keeping the intermediate-tube pressure P\ constant results in higher 
amplitude and higher frequency noise in the test flow. The increase in P l0 results in a lower 
shock speed in the acceleration tube, while the fixed Pi means that the arrival time of the 
reflected expansion will be essentially constant. Thus, at high values of P 10 , the reflected 
expansion is expected to arrive very soon after the shock and introduce high-frequency dis- 
turbances to the test flow. 

To avoid the large-amplitude noise, it should be relatively simple to select operating 
conditions or facility configurations which delay the arrival of the reflected expansion. For 
example, the intermediate tube could be lengthened while fixing the initial gas pressures. 

It is not immediately obvious how to avoid the small-amplitude disturbances observed 
before the arrival of the unsteady expansion. As can be seen in Fig. 23 these disturbances 
do not appear in the test flow for all operating conditions. The acoustic theory of Pauli 
provides some guide to the correct selection of good operating conditions. The basic results 
of that theory are that broad-band noise is focussed to a particular frequency across strong 
expansions and that noise crossing the driver-gas/test-gas interface is attenuated if the speed 
of sound in the test gas is much higher than the speed of sound in the driver gas. Both of 
these effects were observed in the simulations. The stronger driver-gas expansions in cases 
1, 3 and 4 produced significant disturbances in the expanded driver gas while the jumps in 
sound speed across the driver- gas/test-gas interfaces attenuated the noise admitted to the 
expanding test gas. 
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P i 

(kPa) 

Pi 

{kg r/m 3 ) 

Pi 0 

(Pa) 

P 10 

(kg/m 3 ) 

la 

3.45 

5.537 x 10- 3 

16 

2.568 x lO" 5 

ib 

3.45 

5.537 x 10" 3 

16 

2.568 x 10' 5 

2 

20.62 

3.309 x 10~ 2 

16 

2.568 x 10" 5 

3 

1.74 

2.792 x 10" 3 

16 

2.568 x 10" 5 

4 

3.45 

5.537 x 10' 3 

52.6 

8.442 x 10" 5 


Table 1: Initial States. 


Case 

CFL 

cpu — time 
( hrs ) 

memory 
(M words) 

t final 

(ms) 

grid 

la 

0.5 

12.5 

3.2 

5.4 

1604 x 20 

ib 

0.5-0. 6 

42 

6.8 

5.4 

2406 x 30 

2 

0.5-0. 6 

41.6 

6.8 

5.3 

2406 x 30 

3 

0.5 

49 

6.8 

5.4 

2406 x 30 

4 

0.5 

6.5 

3.2 

5.4 

1604 x 20 


Table 2: Computational resources. 


Case 

U.,i 

(m/s) 

U., l0 

(m/s) 

P 5 

(kP a) 

Pi 5 , pitot 

(kP a) 

experimental 
run numbers 

la, lb 

4111 

6914 

1.00 

58.8 

14,134-142 

2 

3375 

6625 

1.01 

122.9 

31 

3 

4203 

6840 

- 

- 

24 

4 

4106 

6100 

2613 

74.7 
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Table 3: Experimentally measured shock speeds and pressures from Shinn and Miller (1978). 


Case 

U, A 

(m/s) 

U s , 10 

(m/s) 

Pb 

(kPa) 

P 5, pitot 

(kPa) 

ib 

3745, 3374 

6820, 5968 

0.867 

51.7 

2 

3108 

6160 

1.067 

127.5 

3 

3711 

5645 

0.713 

33.0 

4 

3490 

5570 

0.733 

72.0 


Table 4: Computed shock speeds and pressures. 


23 












Case 

t sample 

(ms) 

a 2 

(m/s) 

as 

(m/s) 

U 2 
(m/s) 

u 5 

(m/s) 

a 3 
(m/s) 

02/03 

lb 

4.4 

1585 

903 

3351 

6216 

289 

5.49 

2 

4.6 

1594 

616 

2721 

6318 

345 

4.62 

3 

4.4 

1596 

1037 

3567 

6068 

251 

6.37 

4 

4.5 

1620 

1122 

3345 

5443 

281 

5.75 


Table 5: Flow 


properties across the unsteady expansion of the test gas. 


Case 

t rupture 

(ms) 

A U 
(ms) 

A t propagate 

(ms) 

t arrival 

(ms) 

lb 

2.08 

1.00 

2.96 

5.04 

2 

2.40 

1.15 

3.44 

5.84 

3 

2.00 

0.95 

2.82 

4.82 

4 

2.08 

0.95 

2.95 

5.03 


Table 6: Event times associated with the reflection of the unsteady expansion. 
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Primary diaphragm Secondary 

(single or double) diaphragm 

Unsteady 



Figure 1: Conceptual wave diagram for the operation of the NASA Langley expansion tube. 
The numbers refer to the states of the gases as defined by Trimpi (1962). 
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Figure 2: Ideal Pitot pressure history at the exit plane of the acceleration tube. 



Figure 3: Expansion tube geometry as modelled in the simulations. Dimensions are in 
metres. 




x, m 


Figure 4: Wave diagram for case 1. O = shock; A = driver-gas/test-gas interface; + = 
test-gas/acceleration-gas interface; long dash = downstream tail of the unsteady expansion; 
= upstream head of the unsteady expansion; short dash = reflected expansion. 
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°i = 3.45 kPa , P\o = 16 Pa 


primary diaphragm 

t = 0.5 ms 


contact- 


shock 


t = 1.0 ms 


x, m 


1.5 ms 


t = 2.0 ms 


Figure 5: Density contours showing the shock compression of the test gas for case 1. The 
flow states are labelled as per Fig. 1. 










Pi = 3.45 kPa , Pio = 16 Pa 
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Figure 7: Flow properties along the axis during shock compression (case 1, t = 1.0 ms). The 
solid line denotes the ideal one-dimensional solution. 
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Figure 8: Density contours showing the unsteady expansion of the test gas for case 1. 
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Figure 10: Close-up of the flow near the shock and contact surface at t — 4.4 ms, case 1: (a) 
density contours; (b) velocity in a frame of reference stationary with respect to the shock 
and contact. 





6 Pa 
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Figure 12: Comparison of wall pressure and Pitot pressure with experiment (case 1, x 
24 m). The experimental data is taken from Shinn and Miller (1978). 
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x, m 


Figure 13: Wave diagram for case 2. 0 = shock; A = driver-gas/test-gas interface; + = 
test-gas/acceleration-gas interface; long dash = downstream tail of the unsteady expansion; 
o = upstream head of the unsteady expansion; short dash = reflected expansion. 
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Figure 14: Density contours showing the unsteady expansion of the test gas for case 2. 
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Figure 15: Properties along the axis at t = 4.6 ms for case 2: (a) density contours; (b) other 
properties as labelled. 
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Figure 16: Comparison of wall pressure and Pitot pressure with experiment (case 2, x 
24 m). The experimental data is taken from Shinn and Miller (1978). 
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Figure 17: Wave diagram for case 3. 0 = shock; A = driver-gas/test-gas interface; + — 
test-gas/acceleration-gas interface; long dash = downstream tail of the unsteady expansion; 
o = upstream head of the unsteady expansion; short dash — reflected expansion. 
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Figure 20: Wave diagram for case 4. O = shock; A = driver-gas/test-gas interface; + = 
test-gas/acceleration-gas interface; long dash = downstream tail of the unsteady expansion 
_ upstream head of the unsteady expansion; short dash = reflected expansion. 
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Figure 22: Comparison of wall pressure and Pitot pressure with experiment (case 4, x 
24 m). The experimental data is taken from Shinn and Miller (1978). 
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